Environmental effects in the third moment of voltage fluctuations in a tunnel junction 
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We present the first measurements of the third moment of the voltage fluctuations in a conductor. 
This technique can provide new and complementary information on the electronic transport in 
conducting systems. The measurement was performed on non-superconducting tunnel junctions 
as a function of voltage bias, for various temperatures and bandwidths up to 1GHz. The data 
demonstrate the significant effect of the electromagnetic environment of the sample. 
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Transport studies provide a powerful tool for investi- 
gating electronic properties of a conductor. The I{V) 
characteristic (or the differential resistance Rdiff = 
dV/dl) contains partial information on the mechanisms 
responsible for conduction. A much more complete de- 
scription of transport in the steady state, and further in- 
formation on the conduction mechanisms, is given by the 
probability distribution of the current P, which describes 
both dc current /(V) and the fluctuations. Indeed, even 
with a fixed voltage V applied, lit) fluctuates, due to the 
discreteness of the charge carriers, the probabilistic char- 
acter of scattering and the fluctuations of the population 
of energy levels at finite temperature T 0. 

The current fluctuations are characterized by the mo- 
ments of the probability distribution P of order two and 
higher. Experimentally, the average over P is obtained 
by time averaging. Thus the average current is the dc cur- 
rent / = (I(t)), where (.) denotes time average. The sec- 
ond moment (the variance) of P, (SI 2 ), measures the am- 
plitude of the current fluctuations, with SI(t) = I(t) — I. 
The third moment (SI 3 ) (the skewness) measures the 
asymmetry of the fluctuations Q . The existence of the 
third moment is related to the breaking of time reversal 
symmetry by the dc current; at zero bias, 1 = and 
positive and negative current fluctuations are equivalent, 
so (SI 3 ) = 0. An intense theoretical effort has emerged 
recently to calculate the third and higher moments of 
P in various systems 0, 0, IE S H S E3- However, 
until now only the second moment has been measured 
in the many systems studied. Interest in the third and 
higher moments has occurred, first, because its character 
is predicted to differ significantly from that of the second 
moment; in particular, the third moment is insensitive 
to the sample s own Johnson noise for the voltage bias 
case, yet is more sensitive, in a very subtle manner, to 
noise and loading by the environment. Second, measure- 
ments of the higher moments may provide a new tool for 
studying conduction physics, complementary to the sec- 
ond moment. It is the first set of issues that we address. 

In this Letter we report the first measurements of the 
third moment of the voltage fluctuations across a conduc- 



tor, (SV 3 ), where SV(t) = V(t) — V represents the volt- 
age fluctuations around the dc voltage V. Below we relate 
this to (SI 3 ). Our experimental setup is such that the 
sample is current biased at dc and low frequency but the 
electromagnetic environment has an impedance ~ 50 17 
within the detection bandwidth, 10 MHz to 1.2 GHz. We 
have investigated tunnel junctions because they are pre- 
dicted to be the simplest system having asymmetric cur- 
rent fluctuations. However, any kind of good conductor 
can be studied with the techniques we have developed. 
We studied two different samples, at liquid helium, liq- 
uid nitrogen and room temperatures. Our results are in 
agreement with a recent theory that considers the strong 
effect of the electromagnetic environment of the sample 
. Moreover, we show that certain of these environmen- 
tal effects can be dramatically reduced by signal propa- 
gation delays from the sample to the amplifier. This can 
guide future experiments on more exotic samples. 

We present the theoretical overview first for the case of 
voltage bias. In a junction with a low transparency bar- 
rier biased by a dc voltage V , the current noise spectral 
density (related below to the second moment) is given by 
: Sj% = eGVcoth(eV/2k B T) II] (in A 2 /Hz), where e is 
the electron charge and G is the conductance. Only at 
high voltage eV » ksT does this reduce to the Poisson 
result Sj2 = el flU ]. If the barrier transparency is not 
small, the shot noise is Sp = r\el where the Fano fac- 
tor r\ is smaller than 1 1]. For our samples r\ ~ 1. The 
total current fluctuations within a frequency bandwidth 
from h to h is given by (SI 2 ) = 2S P {f 2 ~ fx). The 
spectral density of the third moment of the current fluc- 
tuations in a voltage biased tunnel junction of low trans- 
parency is calculated to be: 5/3 = e 2 GV, independent 
of temperature 0, Q- By considering how the Fourier 
components can combine to give a dc signal, we find that 
(SI 3 ) = 35/3(/ 2 — 2/i) 2 . We have experimentally con- 
firmed this unusual dependence on /i and f 2 - 

We next consider the effects of the sample's electro- 
magnetic environment (contacts, leads, amplifier, etc.); 
the sample is no longer voltage biased. The environ- 
ment emits noise, inducing fluctuations of the voltage 
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across the sample, which in turn modify the probability 
distribution P. Moreover, due to the finite impedance 
of the environment, the noise emitted by the sample it- 
self induces also voltage fluctuations. This self-mixing of 
the noise by the sample is responsible for the Coulomb 
blockade of high resistance samples, where it induces a 
modification of the I(V) characteristics 12] . We con- 
sider the circuit depicted in the inset of Fig. 2, at 
first neglecting time delay along the coaxial cable. The 
noise of the sample of resistance R is modeled by a cur- 
rent generator i. The voltage SV is measured across a 
resistor i?o, which has a current generator in of noise 
spectral density Sp. One has SV = —Rn(i + *o) with 
R D = RRo/{R + R°o) (R in parallel with Rq). The spec- 
tral density of the second moment of the voltage fluctu- 
ations is S V 2 — +R 2 D (Sj2 + Sq). Thus, the only effect 
of the environment on the second moment of the noise 
is to rescale the fluctuations (by R 2 D ) and add a bias- 
independent contribution 0. In contrast, it has been 
recently predicted that the third moment of P is signifi- 
cantly modified by the environment ||: 
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The first term on the right is like that of the second 
moment. The negative sign results from an increasing 
sample current giving a reduced voltage. Our detection 
method is insensitive to S*,3. The environment noise in 
induces voltage fluctuations SV = — Rd^o across the sam- 
ple. These modify the sample's noise Sja (which depends 
on V(t)) as —RoindS^/dV, to first order in SV [l3|. 
This is the origin of the second term. The sample's own 
current fluctuations also modify the sample voltage to 
contribute similarly to the last term of Eq. (1). We give 
below a simple derivation of how to include the effect of 
progation time in the coaxial cable, which dramatically 
affects Sy3. 

Two samples have been studied. Both are tunnel junc- 
tions made of Al/Al oxide/Al, using the double angle 
evaporation technique [l4| . In sample A the bottom and 
top Al films are 50 nm thick. The bottom electrode 
was oxidized for 2 hours in pure 2 at a pressure of 500 
mTorr. The junction area is 15 ^m 2 . In sample B, the 
films are 120 nm and 300 nm thick, oxidation was for 10 
min, and the junction area is 5.6/xm 2 . 

We have measured SV(t) 3 in real time (see Fig. la). 
The third moment of the voltage fluctuations Sy3 is a 
very small quantity, and its measurement requires much 
more care and signal averaging than Sy2 |l5| . The sam- 
ple is dc current biased through a bias tee. The noise 
emitted by the sample is coupled out to an rf ampli- 
fier through a capacitor so only the ac part of the cur- 
rent is amplified. The resistance of the sample is close 
to 50 fl, and thus is well matched to the coaxial cable 
and amplifier. After amplification at room temperature 
the signal is separated into four equal branches, each of 
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FIG. 1: (a) Schematic of the experimental setup, (b) Effect 
of finite bandwidth on the measurement of (SV 3 ). Each point 
corresponds to a different value of the frequencies /1 and fa, 
as indicated in the plot. The data shown here correspond to 
sample B at T = 77 K. 

which carries a signal proportionnal to SV(t). A mixer 
multiplies two of the branches, giving SV 2 (t); a second 
mixer multiplies this result with another branch. The 
output of this second mixer, SV 3 (t), is then low pass fil- 
tered, to give a signal which we refer to as D. Ideally 
D is simply proportional to <SV 3 j where the constant of 
proportionality depends on mixer gains and frequency 
bandwidth. The last branch is connected to a square- 
law crystal detector, which produces a voltage X pro- 
portional to the the rf power it receives: the noise of the 
sample {SV 2 ) plus the noise of the amplifiers. The dc 
current / through the sample is swept slowly. We record 
D(I) and X(I); these are averaged numerically. This de- 
tection scheme has the advantage of the large bandwidth 
it provides (~ 1 GHz), which is crucial for the measure- 
ment. Due to the imperfections of the mixers, D contains 
some contribution of Sya: D — a^Sys + ot^Sy*. We de- 
duce Syz = (D(I) — D(—I))/(2a 3 ), since Sy2 does not 
depend on the sign of I. 

In order to determine the magnitude and sign of (SV 3 ) 
we measured the signal D when the sample is replaced by 
a programmable function generator. The output of the 
generator consist of a pseudo-random sequence of voltage 
steps (at a rate of 10 9 samples per second, lGS/s), the 
statistics of which can be specified. We measured the 
statistics with a 20GS/s oscilloscope. 

Sample A was measured at T — 4.2K. Its total resis- 
tance (tunnel junction + contacts) is 62.6 f2. The resis- 
tance Ra of the junction is extracted from the fit of Sy2 
as a function of eV/ksT, with V the voltage drop across 
the junction. We find Ra = 49.6 fl. Rdiff is voltage 
independent to within 1%. The gain of the amplification 
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FIG. 2: Measurement of S v z{eV/k B T) (solid lines). The 
dashed lines corresponds to the best fit with Eq. (f). The 
dotted lines in the top plot correspond to the different contri- 
butions to Sya (see text). Inset of the middle plot: schematics 
of the equivalent circuit used for the theoretical model. 

chain has been calibrated by replacing the sample by a 
macroscopic 50 VL resistor whose temperature was varied. 
We find rj = 1 with a precision of a few percent for both 
samples. Syz (eV/ftBT) for \V\ < 10 mV is shown in Fig. 
2 (top); these data were averaged for 12 days. 

Sample B was measured at T = 4.2 K, 77 K and 290 
K. The resistance of the junction Rb = 86f2 is almost 
temperature independent. The contribution of the con- 
tacts is ~ 1 57. In Fig. 2 (middle and bottom panels) the 
averaging time for each trace was 16 hours. 

A powerful check that D really measures Sy3 is given 
by the effect of finite bandwidth. Each Fourier compo- 
nent v(f) of SV(t) is amplified by a frequency dependent 
gain g(f) of the setup, such that g(f) = 1 if /1 < |/| < /2 
and g(f) = otherwise. Here /1 and / 2 are positive fre- 
quencies whereas / has arbitrary sign. The measurement 
gives: <^ 3 ) = / (v(f)v(f)v(-f - f')) g(f)g(J')g(-f- 
/') = 35y 3 (/ 2 - 2/i) 2 if f 2 > 2/i, and otherwise. A 
similar argument gives the well known result (SV 2 ) — 
2Sy2 (/ 2 — The scaling of Sy3 with /1 and / 2 is char- 
acteristic of the measurement of a third order moment. 
We do not know any experimental artifact that has such 



behavior [16j. /1 and f% are varied by inserting filters 
before the splitter. As can be seen in Fig. lb, our mea- 
surement follows the dependence on (/ 2 — 2 fi ) 2 , which 
cannot be cast into a function of (/ 2 — /1) |l7l |. 

To analyze our results, consider again the circuit in 
the inset of Fig. 2, a simplified equivalent of our 
setup. i?o ~ 50 n is the input impedance of the ampli- 
fier, which is connected to the sample through a coax- 
ial cable of impedance Rq (i.e., matched to the am- 
plifier) . The sample's voltage reflection coefficient is 
T = (R — R )/(R + R ). In the analysis we present 
next we neglect the influence of the contact resistance 
and impedance mismatch of the amplifier, but we have 
included it when computing the theory to compare to 
the data. The voltage SV(t) measured by the amplifier 
at time t arises from three contributions: i) the noise 
emitted by the amplifier at time t: Roio(t)/2 ; half of «o 
enters the cable, ii) the noise emitted by the sample (at 
time t — At, where At is the propagation delay along the 
cable) that couples into the cable: (1 — T)Ri(t — At)/2 
; iii) the noise emitted by the amplifier at time t — 2At 
that is reflected by the sample: TRoio(t — 2At)/2 ; thus, 

SV(t) = [i (t) + Ti (t - 2At)} - |(1 - T)i(t - At) 

(2) 

For At = 0, Eq. (2) reduces to SV = -R D (i + «o) with 
R D = RR /{R + Rq). Thus, (<5^ 3 ) = -i?^(<z 3 > + 
3(z 2 « ) + 3(?«g) + (i 3 ,)) for A< = 0. In this equation 
the term (i 2 io) leads to the second term on the right of 
Eq. (1). The term (z 3 ) yields the first term of Eq. (1), 
and, due to the sample noise modulating its own voltage, 
the third term of Eq. (1) as well. The terms (zg) and 
(Uq) are zero. The result for At = corresponds to Eq. 
(1), which is a particular case of Eq. (12b) of Ref. 9. 

The finite propagation time does affect the correla- 
tor (i 2 io). The term Sj2 in Eq. (1) has to be replaced 
by + S io( t )i (*-2Ai) 

)/(l + T), where S io(t) 
is the spectral density corresponding to the correlator 
(io(t)io(t — 2 At)). For long enough At this term van- 
ishes, since io emitted at times t and t — 2At are un- 
corrected. Thus, the effect of the propagation time is 
to renormalize the noise temperature of the environment 
To = R S l 2/{2k B ) into T * = T T/(1 + T). 

We now check whether Eq. (1), with Sj3 = e 2 I and 
modified as above to account for finite propagation time, 
can explain our data. The unknown parameters are the 
resistance Rq and the effective environment noise temper- 
ature T * . We checked that the impedance of the samples 
was frequency independent up to 1.2 GHz within 5%. 
Fig. 2 shows the best fits to the theory, Eq. (1), for all 
our data. The four curves lead to Rq — 42 il, in agree- 
ment with the fact that the electromagnetic environment 
(amplifier, bias tee, coaxial cable, sample holder) was 
identical for the two samples. We have measured the 
impedance Z env seen by the sample. Due to impedance 
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FIG. 3: Measurement of S v s {eV/k B T) for sample A at T=4.2 
K (solid lines). AO: no ac excitation (same as Fig. 2, top 
panel). Al: with an ac excitation at frequency £l/2ir such that 
cos 2Q,At — +1; A2: cos 2£lAt = — 1; A3: no ac excitation but 
a 63 Q resistor in parallel with the sample. The dashed lines 
corresponds to fits with Eq. (1). 

mismatch between the amplifier and the cable, there are 
standing waves along the cable. This causes Z env to be 
complex with a phase that varies with frequency. We 
measured that the modulus \Z env \ varies between 30 f2 
and 70 Q within the detection bandwidth, in reasonable 
agreement with Rq = 42 extracted from the fits. 

We have measured directly the noise emitted by the 
room temperature amplifier; we find To ~ 100 K. The 
cable of length ~ 2m corresponds to At long for the 
bandwidth we used. As a consequence, the relevant noise 
temperature to be used to explain the data is T * . For 
sample A, T = 0.11; including the contact resistance and 
cable attenuation one expects Tq = 5 K ; for sample B, 
T = 0.26 and one expects T * = 21 K. A much shorter 
cable was used for T = 290 K, and the reduction of T 
is not complete. These numbers are in reasonably good 
agreement with the values of T * deduced from the fits 
(see fig. 2), and certainly agree with the trend seen for 
the two samples. Clearly T Q * <C Tq for the long cable. 

In order to demonstrate more explicitly the influence 
of the environment on SV 3 we have conducted the follow- 
ing additional measurements. First, by adding a signal 
A sin fit to io (with fl within the detection bandwidth) 
we have been able to modify T * without changing Rq, as 



shown on Fig. 3. The term S 1 ,, 



i (t)i (t- 



-2 At) oscillates like 



cos 2D, At, and thus one can enhance (curve Al as com- 
pared to AO in Fig. 3) or decrease Tq, and even make it 
negative (Fig. 3, A2). The curves A0-A2 are all parallel 
at high voltage, as expected, since the impedance of the 
environment remains unchanged; Rq = 42 fl is the same 
for the fit of the three curves. Second, by adding a 63 f2 
resistor in parallel with the sample (without the ac ex- 
citation) we have been able to change the resistance of 
the environment Rq, and thus the high voltage slope of 



5^3. The fit of curve A3 gives Rq — 24.8 fl, in excellent 
agreement with the expected value of 24.7 fl (63 fl in 
parallel with 42 fl). Since this makes the reflection co- 
efficient negative (r = —0.22), the presence of the extra 
resistor also reverses the sign of Tq , as expected. 

Our data are consistent with a third moment of current 
fluctuations Sj3 being independent of T between 4K and 
300K when the sample is voltage biased, as predicted for 
a tunnel junction. The effect of the environment, through 
its noise and impedance, is clearly demonstrated. This is 
of prime importance for designing future measurements 
on samples with unknown third moment. 
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various mixers, and the asymmetry of S2 vs. current. 
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(8V 3 (I)) measurement. An asymmetric external noise 
would just add just a constant to D, and thus would not 
affect our measurement, even though this constant might 
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Both Sj3 and S V 3 have the same dependence on band- 
width, if the environment impedance is constant. 



